


f 


ns Cain 






Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1999-03 


An infrared model of R/V POINT SUR using 
EOPACE data 


Tan, Chee Yong 


Monterey, California ; Naval Postgraduate School 


http://ndl.handle.net/10945/13653 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 








NAVAL POSTGRADUATE SCHOOL 


t 
Dtic QUA j 


Monterey, California 





THESIS 













AN INFRARED MODEL OF 
R/V POINT SUR USING EOPACE DATA 


by 
Tan Chee Yong 


March 1999 


Thesis Advisor: Alfred W. Cooper 






Approved for public release; distribution is unlimited. 


ji'¥ iNSPHCTE D2 


120 80706661 








Form Approved 


REPORT DOCUMENTATION PAGE 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instruction, searching 
existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this 
burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington headquarters Services, 
Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management 


and Budget, Paperwork Reduction Project (0704-0188) Washington DC 20503. 
2. REPORT DATE 3. REPORT TYPE. AND DATES COVERED 
March 1999 Master’s Thesis | 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 
An Infrared Model of R/V Point Sur using EOPACE data 

6. AUTHOR(S) | 

Tan Chee Yong 


7, PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING 
Naval Postgraduate School ORGANIZATION REPORT 
Monterey, CA 93943-5000 NUMBER 



















1. AGENCY USE ONLY (Leave blank) 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING / 
MONITORING AGENCY 
REPORT NUMBER 


11. SUPPLEMENTARY NOTES 


The views expressed in this thesis are those of the author and do not reflect the official policy or position of the Department of 
Defense or the U.S. Government. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited. 
13.. ABSTRACT (maximum 200 words) 

Infrared polarization techniques to improve the target-background contrast are the subject of much study lately. Polarized 
infrared images of the research vessel, R/V POINT SUR were taken extensively during the Electro-Optical Propagation 
Assessment in Coastal Environment (EOPACE) operational period in Mar-Apr 1996 at Point Loma, San Diego. Contrast 
improvement analysis of these images requires an infrared model of the research vessel. 

This thesis models the apparent infrared signature of the R/V POINT SUR using ship skin temperature records from the 
EOPACE measurements. Using a simple geometric model and adopting a pixel-by-pixel approach, the emitted radiation for each 
pixel surface is calculated from the Planck radiation law. Radiation reflected off the pixel surface is estimated using the polarized 
SeaRad atmospheric propagation code. The total radiance is compensated for atmospheric transmittance with the Navy Aerosol 
Model to arrive at the apparent radiance at the sensor. The apparent radiance is translated into apparent temperature for comparison 
with the recorded images to validate the accuracy of the model. There is moderate agreement between the model and recorded 
images with 45-50% of the pixels falling within 1.75K of the measured apparent temperatures. 


14. SUBJECT TERMS 15. NUMBER OF 


Infrared, polarization, R/V POINT SUR, EOPACE, SeaRad, Navy Aerosol Model Sica 
106 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION OF 19. SECURITY CLASSIFI- CATION 20. LIMITATION 
OF REPORT THIS PAGE OF ABSTRACT OF ABSTRACT 


Unclassified Unclassified Unclassified UL 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 





il 











Approved for public release; distribution is unlimited 


AN INFRARED MODEL OF 
R/V POINT SUR USING EOPACE DATA 


Tan Chee Yong 
Major, Republic of Singapore Navy 
B.A., University of Cambridge, 1990 


Submitted in partial fulfillment of the 
"requirements for the degree of 
MASTER OF SCIENCE IN PHYSICS 


from the 


NAVAL POSTGRADUATE SCHOOL 
March 1999 


Author: 
Tan'Chee Yong 


Alfred W. Cooper, Thesis Advisor 






Ronald J. Pieper, Second Reader 


| : _ _ 


Department of Physics 


lil 





1V 











ABSTRACT 


Infrared polarization techniques to improve the target-background contrast are the 
subject of much study lately. Polarized infrared images of the research vessel, R/V 
POINT SUR were taken extensively during the Electro-Optical Propagation Assessment 
in Coastal Environment (EOPACE) operational period in Mar-Apr 1996 at Point Loma, 
San Diego. Contrast improvement analysis of these images requires an infrared model of 


the research vessel. 


This thesis models the apparent infrared signature of the R/V POINT SUR using 
ship skin temperature records from the EOPACE measurements. Using a simple 
geometric model and adopting a pixel-by-pixel approach, the emitted radiation for each 
pixel surface is calculated from the Planck radiation law. Radiation reflected off the pixel 
surface is estimated using the polarized SeaRad atmospheric propagation code. The total 
radiance is compensated for atmospheric transmittance with the Navy Aerosol Model to 
arrive at the apparent radiance at the sensor. The apparent radiance is translated into 
apparent temperature for comparison with the recorded images to validate the accuracy of 
the model. There is moderate agreement between the model and recorded images with 


45-50% of the pixels falling within 1.75K of the measured apparent temperatures. 
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I. INTRODUCTION 


Every object is a source of infrared radiation with no exception even to the 
atmosphere through which the radiation propagates. Therefore, the design of infrared 
detection and tracking systems requires the ability to discriminate the object of interest 
from a highly cluttered environment. Various filtering techniques such as spectral and 
spatial filtering are typically used to achieve a detectable signal difference between the 
target and the background. More recently, polarization techniques to improve the target to 
background contrast in the infrared wavelengths are being seriously studied. Preliminary 
studies done at the Naval Postgraduate School using the Marine Aerosol Properties and 
Thermal Imager Performance (MAPTIP) trial data showed that the ship and the sky 
background exhibit no significant polarization features while the sea exhibits a 
considerable degree of vertical polarization at near grazing angles. [Ref.1 and Ref.2] The 
studies also showed that horizontal polarization filtering gives a 10 to 20% improvement 
for ship/sea contrast and up to 15% improvement for sea/sky contrast at the horizon. The 
analysis, however, was based on apparent radiance, that is, radiance measured at the 
sensor. It does not show how the relative contribution to the degree of polarization from 
the various sources and atmosphere add up to give the observed polarization features. A 
more conclusive study would need to account for the polarization features based on 
actual radiance derived from radiation theory. Such a study requires good environmental 


data and a reliable infrared model of the ship. 


Electro-optical Propagation Assessment in Coastal Environment (EOPACE) 
organized by NCCOSC-NRaD (now, SpaWarSysCen.San Diego) is a 5-year mullti- 
national effort to improve the performance assessment for electro-optical systems 
operating in a coastal environment. This is because the coastal environment may differ 
significantly from open ocean conditions, and existing propagation codes such as 
LOWTRAN/MODTRAN may not adequately account for the variations. While the main 


effort of the EOPACE measurement series is to fully characterize the coastal 





environment, it also provided the opportunity to study the performance of polarization 


techniques. During one of the EOPACE operational periods in Mar-April 1996 at Point 


Loma, San Diego, extensive polarized infrared images of the research vessel, R/V POINT 


SUR were taken. With the comprehensive environmental database gathered during that 
period, it is possible to perform actual radiance calculations for the ship and the marine 


background. 


This thesis models the apparent infrared signature of the research vessel R/V 
POINT SUR using ship skin temperature records from the EOPACE measurements. The 
research vessel is represented by a geometric model consisting of 23 vertical planes. Each 
plane is associated with a temperature reading from one of the 14 thermistors located 
around the R/V POINT SUR during the imaging period. To facilitate subsequent 
comparison with the recorded images, each plane is broken into pixel size surfaces. For 
each pixel size surface, the emitted radiation from the surface is calculated from the 
Planck radiation law. The radiation reflected off the pixel surface into the sensor is 
estimated from the environmental data using SeaRad [Ref.3]. SeaRad is a background 
add-on to the LOWTRAN 7 /MODTRAN atmospheric codes that calculates the sea 
radiance based on a Cox and Munk [Ref.4] wind driven capillary wave slope distribution. 
In the polarized SeaRad version, the total sea radiance is broken down into its vertically 
and horizontally polarized components. This feature allows us to derive the apparent 
radiance separately for the horizontal and vertical polarizations. For atmospheric 
compensation, the Navy Aerosol Model (NAM) developed specially for marine boundary 
layer propagation is used. Finally, the model is reconstructed into a 2-dimensional array 
with the elements giving the apparent radiance values for each pixel. The array is then 
compared with the ship pixels of the image arrays for each polarization to validate its 
accuracy. A reliable infrared model of the research vessel will offer many more options 
for future analysis, such as contrast improvement analysis, range improvement analysis 


and sensitivity analysis. 








Il. MODELING APPROACH 


A. SURFACE RADIOSITY EQUATION 


Modeling the infrared signature of a target or background is a complex process 


because every object is a source of both emission and reflection radiation. Hence, all 


infrared models are concerned with attempting to solve the surface radiosity equation 


(Ref. 5]: 


L,(6,,,;4) = €,(6,.9,34)B, (7,34) + [6,639,653 4)L,6,.654)dQ, (2.1) 


where L, = 


spectral radiance of the surface S, 

zenith angle of the reflected radiation towards the sensor, 
azimuth angle of the reflected radiation towards the sensor, 
zenith angle of the incident radiation, 

azimuth angle of incident radiation, 

wavelength of the radiation, 

emissivity of the surface, 

Planck 's function, 

surface temperature, 

Bi-directional Reflectance Distribution Function (BRDF) of the surface, 
incident spectral radiance, and 


incident projected solid angle 


This equation essentially states that the spectral radiance of a surface element, S, 


is the sum of two components; the emitted radiation and the reflected radiation. The 


emitted radiation is the product of the Planck function at the temperature of the surface 


and the directional spectral emissivity of the surface. The reflected radiation is the 


integral over all the incident directions of the product of the incident spectral radiance for 
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an elemental projected solid angle of the source and the surface'’s BRDF. 





Sensor 


i 


Figure 2.1. Diagram of geometric relationships [from Ref. 5] 


The surface radiosity equation may seem to be a relatively simple equation to 
solve. However, this is not the case. For instance, there are inter-body reflections. Inter- 
body reflections involve radiation that is reflected to and fro between two surfaces. 
Hence, part of the incident radiance onto a surface may be radiance from the surface itself 
reflected off a second surface back to itself, and that portion of the incident radiance in 
unknown. Furthermore, performing the integral for reflected radiance is not easy and 
usually requires numerical methods. Given the limited time and resources available for 
this thesis, a more practical approach is to simplify the equation to a more useable form 


with a few assumptions. 


Firstly, we see that inter-body reflections are at least twice reflected. In most 
cases, the reflectivity is low and hence, inter-body reflections can be reasonably 
neglected. Secondly, by assuming that the surface is perfectly diffuse, we can use the 


Lambertian approximation to eliminate the angular dependence of the emissivity and 
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reflectivity. Thirdly, assuming that the object behaves like a graybody, and making use of 


Kirchoff's law, we eliminate spectral dependence of the emissivity and reflectivity. 


Hence, we obtain a simpler equation : 


L, (4) = ¢,B,(Z,;4)+(1—e,) (L,@,.4,34) “ (2.2) 
va 


We are still left with the task of performing the integral for reflected radiance, 
which is an angular average of the incident environmental radiance . For each surface 
element orientation, this integral average will be different and will have to be taken into 
account for a collection of elemental surfaces that may represent an object. A further 


simplification is needed and this will be discussed in the later section of this chapter. 
1. Geometric Representation 


To apply the surface radiosity equation to the problem , three basic issues would 
have to be settled first. The first issue is the geometric representation of the object of 
interest. The object of interest is the 135-feet research vessel, R/V POINT SUR, owned 
by the National Science Foundation and operated by the Central California 
Oceanographic Cooperative at Moss Landing Marine Laboratories, Monterey. During the 
recording period, it was used as the cooperative target for infrared imaging. The general 


data of the ship is given in Appendix A. 


The simplest and most convenient way to represent any object is to use a 
collection of flat facets wherever possible. In fact, it is true to say that this is done in 
many infrared models. As observed from the infrared images of the ship at ranges of 
about one kilometer and more, intricate details of the ship are not readily apparent. 
Furthermore, given the relatively small look down angle at the sensor, all horizontal! 
planes have a small viewable area, and do not contribute significantly to the emitted 


radiance. However, it should be noted that at grazing angles, radiance reflected off 


5 

















horizontal planes may be significant. Nevertheless, we will neglect all horizontal planes 
and construct a simple scaled model of the ship using a collection of vertical planes. The 
geometric model is shown in Figure 22a and 2.2b below. The mast structure of the ship 
is not included as it is barely visible in most of the thermal images in this study. As 
mentioned above, the details at the bridge and 01 deck do not show up in the thermal 
images and they are modeled as simple cuboids here. The bow section is approximated as 
2 regular vertical planes (planes 22 and 23). This makes the geometry analysis for the 
bow section much simpler. There are a total of 23 vertical planes, numbered as shown in . 


the figures. The dimensions of each plane are given in Appendix A as well. 


Plane 20 -T15 
Plane 19 -T15 


Sil 
Pe ol 


Plane 21 -T14 


Plane 9 -T8 
Plane 5 -T14 


Plane 22 - T1] | 
Plane | -T10 Plane 2 -T4 Plane 3 -T2 


Figure 2.2a. Port view of the geometric model. and the plane-thermistor matching of 
R/V POINT SUR. The index number such as "T11" refers to the thermistor number as 
listed in Appendix B. 








Plane 17 - T6 
Plane 16 - T6 | Plane 10 - T15 


Plane 12 - T12 Plane 8-T15 


Plane 18 - T5 
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Plane 7 - T16 
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Plane 11 -T 4 Plane 15 - T7 
| Plane 6 - T7 


Figure 2.2b. Starboard view of the geometric model and the plane-thermistor matching of 
R/V POINT SUR 


2. Surface Temperature 


The second issue is to decide how these planes obtain their temperatures. This 
will allow us to calculate the emission term in the surface radiosity equation. One would 
have to consider the sources of heat generation on the ship, such as engines, running 
machineries, electrical switchboards, etc., the various heat transfer processes and also 
interaction with the environment to arrive at a "first principles" temperature distribution 
for the ship surfaces. This is by no means a simple task and is done in more advanced 
models such as Target Contrast Model 2 (TCM2) [Ref.6] adapted for use in the EOTDA 


FLIR Tactical Decision Aid Program. For our purposes, we require a simpler model. 
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For the period of the EOPACE imaging measurements, 14 thermistors were 
attached at selected locations around the R/V POINT SUR. The exact locations of the 
thermistors are listed in Appendix B. The thermistor readings were recorded at 20-second 
intervals throughout the whole imaging period. Hence, we can associate each plane of our 


geometric model with a temperature from one of the 14 thermistor readings. 


In associating the planes to the temperature reading of a thermistor, it is also 
helpful to compare the apparent temperature distribution from the thermal images to the 
thermistor readings. Figure 2.3 below shows a sample thermal image. It can be seen from 
the image that there is a hot section that extends from the stack down to the hull on the 
port side of the ship. As a result, an additional plane was inserted (plane 2) to better 
reflect the temperature distribution of the hull on the port side. (See figure 2.2a) The 
matching of each plane to the temperature reading of a thermistor is indicated beside the 


plane number in Figure 2.2a and 2.2b above. 





Figure 2.3 A typical portside view thermal image from file "Base01.pte". A hot 
section extending from the stack down to the hull of the ship is visible 


3. Radiation from the Environment 


There are two contributors to the reflected radiation term in the surface radiosity 
equation. Since we have ignored inter-body reflections, the environment is the sole 


source of reflected radiation. 








Modeling of the environment is yet in itself another complex process. SeaRad 
with the Navy Aerosol Model is particularly suited for this case since it calculates sea 
radiance, reflected sky radiance, solar radiance, atmospheric path radiance, scattering and 
transmittance over the marine boundary layer atmosphere. However, like LOWTRAN 7 
and MODTRAN, each input case is restricted to a point to point calculation. What this 
means is that the solid angle subtended by the source area has to be reasonably small for 
that calculation to be valid. Therefore, the correct approach is to consider the reflection 
off an elemental surface for each plane. To facilitate the subsequent comparison with the 
recorded images, the size of each elemental surface is made to coincide with the pixel 


size of the image. 


Going back to our simplified equation 2.2, we are still left with the task of 
performing the integral for reflected radiance over the hemisphere of the elemental 
surface exposed to the environment. For each elemental surface, this inte eral would 
require approximately 50,000 input cases to SeaRad. Given the dimensions of R/V 
POINT SUR and the imaging ranges, it is estimated that each frame will have an average. 
of 250 ship pixels. Thus, a huge amount of computational resources is needed for proper 


modeling. A further simplification is required here. 


Consider a pixel surface, which is at a height, h, above the sea level. The range to 
the horizon for the given height h of the pixel is V(Rh) where R is the radius of the earth. 
The angle of elevation of the horizon at the pixel is then given by tan’! V(R/h). (See figure 
2.4 below) Since h is very small compared to R, the value of \(R/h) is very large. This 
gives a to be nearly 90 degrees. Therefore, each pixel will see the sky and sea as 
approximately two equal half hemispheres each subtending a solid angle of 7 steradians. 
If we assume that the incident sea and sky radiance are constant for their respective half 


hemispheres, then equation 2.2 becomes : 


L 
L,(a)= cB(sA)+~6,) “a4 (2.3) 
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Where Le, is the average sea radiance over the lower half hemisphere and 


Lky 1s the average sky radiance over the upper half hemisphere. 






Elemental 
pixel size 
surface 





Elevation angle,a 


Figure 2.4 Geometry relation of a pixel with the sky and sea. Each pixel will see the 
horizon at nearly 90 degrees and hence, the sea and sky as two rougly equal half 
hemispheres. 7 


For a cloudless sky, the sky radiance may be considered to be the sum of the sky 
background emission and the path radiance. It is reasonable to assume that the sky 
background is isotropic. Therefore, sky background emission is independent of look 


angle. However, the path length and hence, the atmospheric path radiance would increase 
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| with increasing zenith angle towards the horizon. Hence, the most appropriate path to be 
used for SeaRad calculation so as to obtain the average sky radiance should be the path 
that gives the average path radiance. This was found to be the slant path to space at a 


zenith angle of 65.3 degrees. 


For the average sea radiance, we will assume the sea to be a flat calm surface. The 
sea radiance is once again given by the surface radiosity equation 2.1. We can write the 


simplified equation as : 
Lisea = Bsea - t(at)(Bsea 7 Lsky) (2.4) 


where r is the reflectance given by the Fresnel equations , 
a is the angle of incidence or elevation and 


we have assumed here that Lsxy is the average sky radiance. 


We can see that Lse, is also a function of the zenith angle. Since Bea and Ly, are both 
constant, the average Le, is then given by the path for average reflectance. The path that 
would best represent the average sea radiance was found to be the slant path to the sea at 
a zenith angle of 118.1 degree. The details of both approximations are given in Appendix 
G. 


ul 
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Ill. IMPLEMENTATION 


The modeling and analysis programs for this thesis were written in Interactive 
Data Language (IDL). However, SeaRad is a FORTRAN code written to run in an MS- 
DOS environment. Therefore, the implementation stage has to be divided into 3 parts, 
switching from IDL to MS-DOS and then back to the IDL environment. 


In the first part, the programs perform 3 main functions. These are : 


- to extract the relevant meteorological data required for SeaRad, 
- to calculate the path geometry between each pixel and the sensor, and 
- to write the meteorological and geometry data into the required input format for 


SeaRad. 


For the second part, the main task is to run the input files through the SeaRad 
code. After running SeaRad, the tasks of the third part are : 


- to extract the required data from the SeaRad output files 
- to calculate the actual radiance received at the sensor 
- to build a 2 dimensional array of the model that will facilitate subsequent 


comparison and analysis. 


A list of the programs written for each part is given in Appendix D. The details of 


operations and information for each part are discussed in the following sections. 


A. PART 1-EXTRACT METEOROLOGICAL DATA, CALCULATE PATH 
GEOMETRY AND WRITE SEARAD INPUT FILES 


There are 4 main procedures written for this part and they run in the sequence : 
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START.PRO, PLANEIN.PRO, SEAIN.PRO, AVESEAIN.PRO. Each main procedure 
performs a specific function. START.PRO extracts the relevant meteorological data. It 
also provides the basic geometry parameters to the procedure PLANEIN.PRO, which 
calculates the pixel coordinates and pixel-sensor path geometry. After all the required 
data are extracted and calculated, the procedures SEAIN.PRO and AVESEAIN.PRO are 
used to write the data into FORTRAN input format required to run the SeaRad code. The 


details given in the subsections below will more or less follow this sequence of operation. 





1. Data File Format 


The original infrared images of the R/V POINT SUR were recorded using the 
PTRWIN, version 3.16 data acquisition software from CEDIP. They are identified with 
the file extension .PTW. Although PTRWIN has the capability to analyze the images, it is 
inadequate for analysis with large numbers of images. Therefore, they were converted 
into a more efficient format using an IDL procedure written by Pontes [Ref. 7] called 
BASEGEN.PRO. The converted files, identified by the filename "BASE*.PTE", have a 
basic structure that contains an information vector, a time vector that gives the time of 
each frame and is followed by the images in the corresponding order. The detailed 
structure and type of information carried by each "BASE*.PTE" file is given in Appendix 
E. 


Each "BASE*.PTE" carries a total of approximately 75 picture frames, that is, 
roughly 25 frames for each case of horizontally polarized, vertically polarized and 
unpolarized images. Given that the scanner produces 6.25 picture frames per second and 
that each image set of different polarization is taken sequentially, there is a time lag of 
about 8 seconds between the first frame and the last frame of the file. One of the 
objectives in the previous analysis performed by Pontes [Ref. 7] on the EOPACE data 
was to determine whether there is a ship radiance variation with time. It was found in that 
analysis that the variation of radiance within the time scale of 8 seconds is not significant. 


Therefore, only the first frame of each polarization case will be considered in this thesis. 
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2. Geometry of Measurements 


The measurements were made from the NCCOSC Building 15 (32° 39' 55" N, 
117° 14' 32" W, elevation above sea level -57 feet) at the tip of Point Loma. The research 
vessel was imaged at different aspect angles at ranges from 0.5 km up to 5 km. The 
bearing of the ship was also made to vary from about 100 to 260 degrees from the 


observation point. A chart showing the location of the measurement is in Appendix F. 
3. Meteorological Data 


Meteorological data were collected at 2 separate meteorological towers (METOC 
1 and METOC 2) onboard the R/V POINT SUR. The data from METOC 1 and METOC 
2 were recorded in files identified by "EOP*.TXT" and "EOP30SEC.*" respectively, 
where * indicates the day of the year. Each file was started at 0:00 Greenwich Mean 
Time (GMT) of each GMT day and the data were recorded in 30-second intervals. For 
each 30-second interval, 23 parameters were recorded. These are the GMT time of the 
reading, relative wind speed and direction, ship speed and direction, true wind speed and 
direction, air temperature, relative humidity, pressure, sea surface temperature and the 
GPS position of the ship. The detailed format of the METOC files is included as 
Appendix G. 


To extract the meteorological data, the time vector in the .PTE file is first read. 
This is then converted to GMT to open the correct files and to read the nearest line of 
data recorded for that time. Since not all the data are needed, the relevant data are 
extracted, averaged and used to calculate some additional required information. The 
program READMET.PRO does all these, and returns a vector containing the range of 
ship, sea temperature, air mass parameter, current wind speed, 24 hour wind speed, wind 
azimuth relative to observer azimuth, GMT day of the year, GMT and observer azimuth. 


This information for the analyzed files is documented in Appendix H. 
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The US Navy Aerosol Model (NAM) was developed to compute the aerosol 
extinction coefficients for visual and infrared wavelengths in the marine boundary layer. 
Since we are concerned with the propagation of infrared radiation from a ship in a coastal 
environment , it is particularly suited for application in this thesis. However, in order to 
use the Navy aerosol model, part of the information required is the "air mass parameter 
(AMP)", which is included in the MODTRAN propagation code as "ICSTL". This 
parameter is not recorded in the METOC files and cannot be derived from any of the 


parameters collected there. 


The traditional way to determine the air mass parameter is to measure the gas 
concentration of the radioactive Radon-222 in the atmosphere. During the Intensive 
Operational Period (IOP) in Mar-April 1996, the Radon-222 gas concentration was 
measured at San Nicolas Island, Point Mugu and San Diego. This information is available 
to the EOPACE participants at the EOPACE website at http://sunspot.spawar.navy.mil:80 
/343/eopace/sd_iop_ma96/nrad883/bdg323/radon/radonsd.txt". With these 


measurements, the air mass parameter is then simply given by the equation [Ref.9] : 


R 
P= wv1| =| i (3.1) 


where INT truncates to the integer value and Rn is the atmospheric Radon-222 gas 
concentration in picocurries per cm*. It should be noted here that this equation is 
empirical and would have to be modified for application in other geographical regions. A 


brief discussion on the NAM and the air mass parameter is in Appendix I. 
4. Pixel Size 
The infrared scanner used to record the ship images is the dual band serial 


scanning AGA-780 Thermovision Thermal Imaging System. The short-wave (SW) 
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scanner uses an indium antimonide (InSb) detector and covers the 3-5.6 jm spectral 
band. The long-wave (LW) scanner covers the 8-14 um spectral band and has a mercury- 


cadmium-telluride (HgCdTe) detector. 


Electro-magnetic radiation from an object is focused onto the detectors by two 
separate systems of infrared lenses and rotating prisrhs. In each scanner, an infrared 
image frame is produced by 2 prisms, rotating in such a way that they produce four fields 
to make one interlaced frame. This is shown in Figure.3.1 below. Each field has 100 
horizontal scanning lines and only 70 of these are used as active imaging lines. This 


gives a vertical scanning efficiency (Msc) of 70%. 
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Figure 3.1 Interlace fields of the AGA-780 [Ref.9] 


There were 3 different lenses that were available for use with the scanner. These 
are the 3.5°, 7° and 20° field of view (FOV) lenses. Table 3.1 lists the specifications for 


these lenses. 
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Table 3.1. Lens specification from AGA 780 Manual [Ref.9] 





For the whole recording period during EOPACE April 96, only the 3.5° and 7° 
lenses were used. The horizontal pixel resolution is therefore 0.5mrad and 1.1 mrad 
respectively. To obtain the effective vertical pixel resolution, we multiply the geometric 
resolution by the scanning efficiency (ns-) . This gives 0.35mrad and 0.77mrad for the 
3.5° and 7° FOV respectively. 


5. Path Geometry 


The idea behind the path geometry calculation was to first determine the 
coordinates of the center of each pixel (which are called nodes) with respect to an origin 
vertically below the sensor. Thence, the sensor to pixel path is a vector given by the 


difference of the pixel coordinates and sensor coordinates. 


To obtain the 3 dimensional coordinates of each node, we need to consider each 
plane of the geometric model of the ship (Figure 2.2, Chapter IT) one at a time. The 
program PLANE.PRO is first used to select the visible planes for computation. This 
significantly reduces the computation load. Each selected plane is then broken down into 
nodes depending on the size of the elemental surface, which in turn, is a function of the 
range and pixel resolution. The position of each node relative to the center of the ship is 
first calculated. A rotation is then applied to the node positions to orientate the plane to 
the correct aspect angle. Lastly, the origin is shifted from the center of the ship to an 


origin directly below the sensor. 
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One of the main problems at this stage is the fractional edges of the planes. As 
each plane may not contain whole number of pixels, the edges have to be carefully 
handled so that when the neighboring plane is put together, it does not give an incorrect 
radiance value later on. To do this, the pixels at the edge are always rounded up. Then a 
weighting factor based on remaining area at the edge to the pixel area is calculated. For 
those pixels at the edge of a plane, the weighting would be less than one while the 
weighting for the pixels at the edge of the adjacent planes will be one minus the 
weighting of the edge pixels of the first plane. These pixels will share the same node 
positions. Therefore, the output of the routine PLANEIN.PRO contains a weighting 


matrix in addition to the 3 dimensional coordinates of the position of the nodes. 


The reflection path was discussed towards the end of Chapter II. The paths 
representing the average sea and sky are just the slant paths with zenith angles of 118.1 


and 65.3 degrees at the pixels respectively. 


6. SeaRad [Ref. 3] 


SeaRad was developed as an ocean background add-on to the MODTRAN 2 and 
LOWTRAN 7 atmospheric codes. The code used to calculate atmospheric transmittance, 
solar, sky and path radiance is the same code as MODTRAN/LOWTRAN. The sea 
radiance add-on is activated by a "seaswitch". When this switch is set to true, the code 
will take the boundary temperature as the sea surface temperature whenever the path 
intersects the sea surface. It then applies the Fresnel coefficients to the reflection of 
infrared radiation by a seawater capillary wave facet. The result of a single tilted facet is 
then integrated over all possible tilts which follows a Cox and Munk probability 
distribution function. Hence, SeaRad calculation gives a pixel average sea radiance, 


which suited our purpose here. 


In addition, it is worth noting that SeaRad does not take into account the 
correction for shadowing effects proposed by Saunders [Ref. 10]. In essence, Saunders 


reasoned that radiation from certain portions of a wave is blocked by the wave in front. 
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Hence, these portions do not contribute to the sea radiance. The amount of area being 
blocked is dependent on the viewing angle. The more acute the viewing angle, the more 
area of the wave is blocked. Given the viewing geometry of our problem, this effect 


maybe appreciable. 


SeaRad input requires the submission of a series of at least 5 cards. Each card isa 
line of formatted input. As each pixel requires one set of input cards and there is a large 
number of pixels, the routines SEAIN.PRO and AVESEAIN.PRO are used to write all 
the extracted meteorological data and calculated geometry data into the required format 
for each pixel. A sample input file is shown at Appendix J. An explanation of the input 


_ cards used for the input files is also included in Appendix J. 


B. PART 2 - RUN SEARAD 


SeaRad is run in an MS-DOS environment. The DOS commands to run a single 


input file are : 


copy "input filename" tape 5 
SeaRad 


The output of the calculations is written to an output file called OUT. As can be 
seen, SeaRad is not designed to run a large number of input files as is needed here. Since 
each input file differs only in the geometry data, it would be possible to modify the 
SeaRad codes to accept multiple geometry inputs. A simpler alternative, though less 
efficient, is to run SeaRad in batches using MS-DOS batch file commands. A batch file 
called SEA1.BAT was written to run input files in batches. There are 3 types of input 
files in the "CASEI/INPUTS/" directory : the first type is for the pixel-sensor path | 
geometry , the second type for the pixel-sea geometry and the last type for the pixel-sky 
geometry. These files are identified by the filename "*p*c**.std", "*a*b*0.std" and 
"*a*b*1.std" respectively. The output of each file has exactly the same name as the input 


file but they are copied into the "CASE1/OUTPUTS(/" directory. 
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C. PART 3 - EXTRACT SEARAD OUTPUT, CALCULATE RECEIVED 
RADIANCE AND CONSTRUCT 2-D ARRAY 


After running SeaRad, the procedure CONTINUE.PRO reads the SeaRad output 
files to extract the required information. It then calculates the radiance received at the 
sensor and reconstructs the model into a 2 dimensional array. The second section of | 
CONTINUE.PRO compares the resultant array with the recorded image arrays. The 


second section will be discussed in the next chapter. 
1. SeaRad Output 


The format of the SeaRad output files depend on the type of calculation that is 
done. For the three input cases, two types of output files were written. A sample of each 
type of output file is shown in Appendix K. The routines SEAOUT.PRO and 
AVSEAOUT.PRO are used to extract the required information from the output files. The 
information returned by these programs is : 


- pixel to sensor path transmittance ( Tpath ), 

- pixel to sensor path radiance ( Lpatn), 

- total radiance at the pixel for the average sky radiance path (Ls) 

- total radiance at the pixel for the average sea radiance path in horizontally 


polarized (iiss) vertically polarized (Lysea).and randomly polarized components 


(Lusea). 
2. Ship Skin Temperature 


As previously mentioned, the ship skin temperatures of R/V POINT SUR were 
recorded continuously in 20-second intervals throughout the whole measurement period. 
These readings were stored in a file called "SKINDATA.PTE". Each line has 16 


readings corresponding to 2 reference temperature readings and 14 thermistor readings. 
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The detailed format for the skin temperature file is in Appendix L. 
3. Ship Skin Radiance 


Given the ship skin temperature, the spectral radiant exitance from the ship skin 


surface can be calculated using Planck radiation law : 


2ahe’ 
M 5,(A,T) = 2 


Fea) (3.2) 


where T is the surface temperature in K, 
Xd is the wavelength in m 
h is Planck's constant (6.6256 x 10 Js), 
c is the speed of light (3 x 10° m/s), 
k is Boltzman's constant (1.38054 x 10°? J/K) and 


Mpz has dimensions in watts per meter square per wavelength interval in m. 


The radiant exitance for a graybody over x steradians of solid angle is then given 
by the Stefan-Boltzman law, which is : 
j (3.3) 
M oraysoa(T) = € { M,,(A,T)dA = eoT 


where ¢ is the emissivity of the surface and 


o is the Stefan-Boltzman constant (5.6697 x 10°)? W/em’K’). 


However, we are only interested in the in-band radiance for the long-wave band 


(8-12 um). The in-band radiance can be represented as : 


& A2yun 
Lape =~ |” M gp (4,7 aA (3.4) 
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There are two methods to calculate the in-band radiance. One method is to use the 
universal blackbody curve or pre-tabulated integrals over it. The other method is to 


perform the integration numerically. The routine "EMISSION.PRO" performs the 


numerical integration using step size of 0.01 um. 


Since we have already associated each plane with the temperature reading of a 
thermistor, the in-band radiance is calculated for each pixel using the above equation. In 
the calculations, the ship surface is assumed to behave like'a graybody with an emissivity 
previously measured to be 0.95. Since no further information on the optical properties of 
the ship paint is available, another assumption was made that the horizontal and vertical 


polarized components of the in-barid radiance are equal. 
4. Radiance Received at Sensor 


Finally, the radiance received at the sensor is the sum of the ship skin radiance 
compensated for atmospheric transmittance, the path radiance and the reflected sea and 
sky radiance. This is best summarized by the following formulae for the two polarized 


and randomly polarized components : 





| L, 

Licensor ee 3 +L path + (I = ENT yath Le a “2 (3.5a) 
1 L, 

Lisi = 3 +L path + (I i Vat Le | (3 Sb) 
| (1-<) 

Lijsnenp = ery Epa a3 L path By 9) tv path (= + i) . (3.5c) 


It should be noted here that there would be a transmission loss due to the polarizing 
filters. However, it is not necessary to correct for it since the corrections are built into the 


scanner calibration process. 
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5. Reconstructing a 2-D Array 


Remembering that the procedure was done one plane at a time, the remaining task 





is to put all the planes together to give a 2 dimensional image of the model. This process 


is done in 3 steps. 


The first step is to multiply the radiance matrix of each plane by the weighting 
matrix so that the pixels at the edges of adjacent planes will fit together. The second step 
is to correctly position each plane relative to a reference point (the origin at the midpoint 
of the ship) in an empty 72 x 30 matrix. This is done by the procedure PLACE.PRO. An 
additional procedure called BLOCK.PRO is then applied to blank out the areas of overlap 
between planes for that viewing angle. Lastly, the radiance values of the planes are | 
simply added together to form a 2D array of ship radiance received at the scanner. A 
sample image(magnified) of the ship model derived for file "Base01 pte" is shown in 
Figure 3.2 below. The recorded image is the original image from the data files. The 
extracted image is the resultant image when a mask having the same shape as the model 


but containing ones and zeros is applied to the original images. 





Recorded image 





Model Extracted image 


Figure 3.2 Images of R/V POINT SUR for unpolarized case from file Base01.pte. The 
model image shows a higher temperature than the recorded temperature in the extracted 
image. Notice that some of the background pixels are inevitably included in the extracted 
image 
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IV. COMPARISON OF RESULTS 


A. APPROACH FOR COMPARISON 
1. Pixel-by-pixel Comparison 


Since we started out using a pixel by pixel approach to derive a 2-dimensional 
array for received radiance at the sensor, each element in the 2-D array will directly 
correspond to a pixel in the recorded image frame. Each image frame has a sea and sky 
background. Since we are only interested in comparing the ship model with the measured 
ship radiance the background would have to be eliminated. To extract the ship pixels in 
the image frame, a mask is first created from the 70 by 30 model array. This is done by 
setting elements of ship pixels in the array to ones and the background to zeros. Next, the 
coordinates of the right bottom corner of the ship in the image frame are read from user's 
input. This information is used as the reference point for aligning the mask over the 
image frame to extract the ship pixels in the image. Having set all the background 
elements of the image array to zero, we can now compare the two arrays directly by 
simply taking their difference. However, before the difference can be taken, the elements 
in the arrays must have the same units. The model array has units of radiance while the 


image array has units of digital level. A unit conversion is required. 
2. Thermal Level 


The numerical measure of the received and detected radiation used in the infrared 
scanner is called the thermal level. The traditional unit for the thermal level is the 
isotherm unit, which is linearly related to the received photon radiation. Therefore, the 
relationship between the thermal level and the object temperature is non-linear. To 
translate object temperature to its thermal level, a calibration function is used. This 


calibration function is obtained by plotting the thermal level against the known 
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temperature of a blackbody. It has a characteristic equation of a form similar to the 


Planck equation : 


TL = —4A___ cs (4.1) 


where A, B, C are constant determined by curve fitting, 
OS is the offset correction for detector drift and 


T is the temperature in K. 


The calibration function is dependent on the aperture and filter used and a set of 
calibration curves is usually produced. Prior to the measurement period, the scanner with 
polarizing filters was calibrated. The calibration curves are shown in Appendix M. The 
calibration constants for each polarization filter are also given in the appendix. Although 
all the calibration constants are given, it should be noted here that the offset constant for 
detector drift, OS, tends to vary rapidly with changing ambient conditions. Using the 
calibration functions, the thermal level of the ship pixels for each polarization case in the 
images can be easily translated to the equivalent blackbody temperature of the ship 


surfaces. 
3. Digital Level 


It was found that the images recorded by the PTRWIN software were stored in 
digital levels in the original data files. Although the PTR WIN software have the function 
to convert digital levels to isothermal units, the process to convert all the "*.PTE" files 
into thermal levels would require re-storing the original "*.PTW" file into isothermal 
units and re-generating the "*.PTE" files using BASEGEN.PRO. An alternative way is to 
deduce the relationship between the digital level and thermal level from the PTRWIN 
software. It was found that they are related by a simple linear relation which depends on 
the Range (R) and Level (L) settings in the PTRWIN software. Hence, to convert digital 
level (DL) to thermal level (TL), the following empirical formula is used : 
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TL=A*DL+B (4.2) 


where A=0.00122, B=L-2.5 forR=5, 
A=0.00244, B=L-5 forR=10 and 
A=0.00488, B=L-10 forR=20. 


The range and level settings are not available in the information vector of data 
(* PTE) files. They have to be extracted from the original data (*.PTW) files. Appendix 


N documents the original files and their range/level settings. 
3. Calibration Drifts 


The time lag between the first frames of the horizontally and randomly polarized 
images in a single data file is roughly g seconds. Within this time scale, the temperature 
distribution of the ship surfaces has not changed significantly. If we translate the first 
image frame of each polarization in a data file into the equivalent blackbody temperature 
using the calibration curves, we should expect the blackbody temperature arrays for the 
three different polarizations to be almost similar, that is, a pixel in one polarization frame 
should have the same equivalent blackbody temperature as the corresponding pixel in 
another polarization frame. Hence, if we take the difference between the blackbody 
temperature arrays of any two polarizations, the elements of the resultant temperature 
difference matrix would be mostly zeros. In another way, the average value of all the 
elements in the resultant temperature difference matrix should be nearly zero. Table 4.1 
below shows the average values of the elements of the resultant temperature difference 
matrices when the blackbody temperature arrays of the three different polarizations of a 
single data file are compared to each other. Column one gives the average temperature 
difference between horizontal and vertical polarizations: column two , horizontal and 


random polarizations; and column three vertical and random polarizations. 
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Table 4.1 Mean pixel temperature difference between polarization cases of an image 


file 


Average temperature difference (K) 
Base09.pte -0.053289222 -6.0796106 -6.0263214 


We can see that the temperature differences between corresponding pixels of the 





















polarized image and unpolarized image are significant. If this difference were solely due 
to the contrast improvement of the polarizing filters, we should expect opposite signs for 
the second and third columns. Clearly, this is not the case and it is likely that calibration 
drifts are present. Since we do not know exactly what the calibration drifts are, a relative 
calibration is required in order for us to relate the measured ship skin temperatures to the 
recorded thermal levels. The procedure for a relative measurement is described in section 


10 of the AGA manual [Ref.10 ]. The relevant extract is reproduced in Appendix O. 
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To perform the relative calibration, we use the following formula from the AGA 


manual : 
24+, (4.3) 


where I, is the relative calibrated thermal level, 
I, is the calibrated thermal level of a reference point on the object whose 
temperature is known, | 
Al,, is the thermal level difference from the reference level (1;), 
T>. is the atmospheric transmittance between object and scanner and 


€ is the object emissivity. 


Since the values of the elements of the matrix that we derived are radiance values 
received at the scanner aperture, these would have to be converted into equivalent 
blackbody temperatures. The routine RADTOT.PRO converts radiance to equivalent 
blackbody temperature by linear interpolation in three iterations. Since the radiation- 
temperature relation is actually non-linear, a good initial guess of the temperature range 1S 
required to achieve good conversion accuracy. For the temperature range of the ship, the 


conversion accuracy is 0.001 W/m’. 


If we now use the converted blackbody temperature of the pixels for the 
recalibration, then t, and € are simply equal to one in equation 4.3 since these are the 
equivalent blackbody temperatures at the aperture of the scanner. In addition, any pixel 
can be used as the reference point since now they each have a known equivalent 


blackbody temperature. 
4. Unit of Comparison 


The comparison can be done using radiance, temperature or thermal level 
differences. However, thermal level differences do not give a good indication when we 
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compare across polarization cases. For example, a thermal level difference of 1 in the 
horizontal polarization case may be due to-0.3 W/m? difference in radiance while a 
difference of 3 isotherm units in the unpolarized case is just 0.15 W/m’. Hence, doing the 
comparison in temperature units is the most convenient since the recalibrated thermal 


levels of the image pixels are easily translated to temperature. 
B. DISCUSSION ON RESULTS 
1. Statistical Measures 


In an ideal case, the temperature difference between each pixel of the model and 
the recorded image would be zero. More realistically speaking, we would expect a 
distribution of pixels with a mean of close to zero and a small standard deviation. Figure 
4.1 below shows the pixel distributions for a sample file. From the figure, a large number 
of pixels are clustered close to the zero temperature difference mark. However, there are 
a small percentage of pixels that register a big temperature difference. The origins of | 


these outlying pixels are likely due to two limitations. 


The first limitation is inherent in the geometric representation of the ship. As our 
geometric model is a simplified representation of the ship, not all calculated pixels will 
have a corresponding pixel in the image or vice versa. Secondly, the alignment of the 
mask is highly sensitive, that is, down to pixel resolution. At this level of resolution, it is 
difficult to achieve a 'correct' one to one correspondence in the ship pixels. The impact of 
all this is that some pixels will be compared with the background pixels resulting in a big 
temperature difference. It has the effect of driving up the mean and standard deviation 


especially when the number of ship pixels is small. 


In order to obtain a better characterization of the pixel distribution, two additional 
statistical measures were calculated. These are the modal temperature difference and the 


fraction of pixels falling within a specified range from the modal temperature difference. 
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The modal temperature difference has the advantage that it is not affected by the outlying 
pixels. However, it is more susceptible to freak occurrence of peaks in a distribution. 
Therefore, both sets of statistical measures are needed to gain a fair interpretation of the 


results. 
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Figure 4.1 Pixel distribution of difference in equivalent blackbody temperature 
between model and recorded image for file Base10.pte. The small percentage of 
pixels at the outlying areas tend to drive up the mean and standard deviation 
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2. Results of Comparison 


A total of 37 data files were analyzed and compared. The results of each 
comparison are tabulated in Table 4.2. Most of the cases analyzed showed a modal 
temperature difference of less than 1K between model and image. The average 
temperature differences for these cases were also less than 1K. The modal temperature 


difference, in general, is lower than the average temperature difference. 


There are 5 files where the modal temperature difference is much higher than the 
average temperature difference. The modal temperature differences for the 3 polarization 
cases in these files were very different as well. In contrast, the average temperature 
differences for the 3 polarization cases were quite consistent. It is likely that there are 
occurrences of freak peaks in the distribution making the modal temperature difference 


unreliable. 


There is one file that has average and modal temperature differences of more than 
2K for all 3 polarization cases. This is file Base07.pte. The total number of pixels for 
each frame of this file is 32. It is likely that there are too few good pixels to give a 
reasonable comparison. With the exception of file Base07 pte, there is moderate 
agreement between the model and recorded images with 45-50% of the pixels falling 
within 1.75K of the measured apparent temperature. | 
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File no. Number | Polarisation Measure of centre of Measure of spread 
of pixels case distribution 


(base*.pte) 


Modal Mean Fraction of pixels | Fraction of pixels | Standard 
temperature | temperature | within 0.25 K of | within 0.75 K of | deviation 
difference difference modal temp. modal temp. 
difference difference 
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File no. Number | Polarisation Measure of centre of Measure of spread 
of pixels case distribution 


Modal Mean Fraction of pixels | Fraction of pixels | Standard 
temperature | temperature | within 0.25 K of | within 0.75 K of | deviation 
difference | difference modal temp. modal temp. 
difference difference 


(base*.pte) 
15 92 -0.945 -0.209 0.207 0.511 1.917 
-0.164 -0.070 0.228 0.505 2.219 
-0.931 -0.932 0.228 0.467 1.793 
"|" f= -0.945 0.569 0.153 0.378 2.64 
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File no. Number | Polarisation Measure of centre of Measure of spread 
of pixels case distribution 


(base*.pte) 


Modal Mean Fraction of pixels | Fraction of pixels | Standard 
temperature | temperature | within 0.25 K of | within 0.75 K of | deviation 
difference difference modal temp. modal temp. 
difference difference 
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File no. Number | Polarisation Measure of centre of Measure of spread 
of pixels case distribution 


(base*.pte) 


Modal Mean Fraction of pixels Fraction of pixels | Standard 
temperature | temperature | within 0.25 K of | within 0.75 K of | deviation 
difference difference modal modal 
temp.difference. | temp.difference 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The simplified surface radiosity equation is the basis for the modeling. It states 
that for each pixel the radiance from the surface is the sum of the emitted radiance and 
reflected radiance. Since we have made the Lambertian and graybody approximations, 
the emitted radiance accounts for almost 95% of the radiance before atmospheric 
attenuation. From the comparison of apparent temperatures between the model and the 
images, 55% of the pixels have temperature deviations of more than 1.75K from the 
images. For the range of temperatures read by the thermistors, this deviation in radiance 
is of the order of about 3 Whm2-sr. We can see from the SeaRad outputs (Appendix K) 
that the average sea and sky radiance incident on the pixel surface is of the order of 30 
W/m2-sr and 10 W/m7*-sr respectively. This gives the reflected radiance to be in the order 
of 1 W/m*-sr. Therefore, the origin of the 1.75K temperature deviations for 55% of the 
pixels cannot to be due to the reflected radiance. A second possibility for the deviation 1s 
the path radiance due to path length error. The path radiance is of the order of 0.01 W/m’- 
sr per meter increase in range. It would require a range error of more than 300m to 
account for the 3 W/m*-sr deviation. Again, this is unlikely to be the source of deviations 
since the range or elevation error is not expected to be more than 10m. Therefore, we can 
conclude that the mostly likely source of the error is the surface temperature distribution 
used in the model. It does not sufficiently represent the actual ship skin temperature 


distribution. 


One of the problems with modeling the temperature distribution of the ship is that 
there are too few thermistor readings. As a result, each plane is assumed to have a | 
constant temperature associated with the nearest thermistor. Given the limitation on the 
number of thermistors, one option is to perform sensitivity analysis with the model. This 
will help in a better understanding of the actual temperature distribution of the ship. One 
approach is to perform a plane by plane analysis. Additional planes can be added to give 


a better distribution where necessary. Another approach may be to consider one or two 


4] 














dimensional heat flow for each plane so as to arrive at a temperature distribution for the 


pixel surface elements. 


As discussed earlier, the main contributions to the observed apparent radiance are 
the emitted radiance from the ship surfaces and the path radiance. Since we have assumed 
these two contributions to be randomly polarized, the apparent radiance that we derived 
from our ship model would show no significant polarization features. When our model is 
compared to the measured images for the three polarization cases, the match is 
comparable in all three polarizations. The average overall modal temperature differences 
are 0.078K, -0.093K and —0.093K for horizontal, vertical and random polarization 
respectively. The spreads of the distributions for all three polarizations are also 
comparable at 40-45% of the pixel within 0.75K of the mode. This shows that the 


assumptions used in our model are justified. 


For the whole modeling process, the running of SeaRad code is the most time 
consuming. This is because each input file is copied to the standard input file to run the 
SeaRad code and the standard output file is then copied to an output file. As each pixel 
requires one input file, it takes approximately one hour to run a frame consisting of 360 
pixels. It also takes up a lot of memory space to store the input and output files. The 
SeaRad input data requirement can be broken down into two types of data. The first type 
is the meteorological data. This data does not change for one frame and is the same for all 
input files of that frame. The second type is the geometry input. This is the part that | 
changes for each pixel. If it were possible to modify the SeaRad code to accept multiple 
geometry input, then this would significantly reduce the processing time as well as the 


time spent copying files. It will make the SeaRad code a much more practical tool. 
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APPENDIX A. DIMENSIONS OF PLANES IN GEOMETRIC 
MODEL AND GENERAL DATA ON R/V POINT SUR 


A. DIMENSIONS OF PLANES IN GEOMETRIC MODEL 


The geometric representation of the ship is made up of 23 vertical planes. An 
illustration of the ship model is at Figure 2.2 in Chapter 2. The table below summarises 


the dimensions of each plane in the model. 


Table Al. Summary Of Planes In Ship Model 


rs. 







Matched to 






thermistor 
Port hull below stack 
01 deck port side 
Stack port side 


Starboard hull aft section 
12.6 4.5 Starboard hull mid section 
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thermistor 
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APPENDIX B. THERMISTOR LOCATIONS ON R/V POINT SUR 


The thermistor location is described at Table B1. Figure B1 shows the locations. 


Table B1. Location of thermistors 
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Figure B1. Location of thermistors 
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APPENDIX C. AVERAGE SEA AND SKY RADIANCE 
APPROXIMATIONS 


A. AVERAGE SKY RADIANCE PATH 


The sky radiance may be considered to be the sum of the sky background 
emission and the path radiance for a cloudless sky. It is reasonable to assume that the sky 
background is isotropic and therefore, sky background emission is independent of look 
angle. However, path radiance would increase with increasing zenith angle towards the 
horizon since the path length is longest at the horizon. Hence, the path representing the 


average sky radiance is the same path as that for the average path radiance. 
To obtain the average path radiance, we consider the geometry of the atmospheric 


layer and the earth surface in Figure C1 below. SeaRad assumes an atmospheric layer of 


100 km from the earth surface. 


R +100 


Figure C1. Relationship between zenith angle, 0 and angles a and B 


49 





The angles 8, a and B are simply related by the sine rule : 








sinl80-@) sinf sina (C.1) 


Therefore, we find that the path length, P is a function of 0. The average path 
length is then given by: 


7 276 (C.2) 
y= JP@ae 


Performing the integral numerically, we obtain an average path length of 208 km, 


which corresponds to a slant path to space at a zenith angle of 65.3 degrees. 
B. AVERAGE SEA RADIANCE PATH 


Assuming the sea to be a flat calm surface, the sea radiance may be simplified as: 


Lsea = Bsea (T)- p(9; )[Bsea (T)- Lsky] (C.3) 


where p is the reflectance given by the Fresnel equations , 
0; is the angle of incidence or elevation, 
Bsea is the blackbody emission of the sea surface at temperature T and 


we have assumed here that Ly is the average sky radiance. 
We can see that Lea 1s also a function of the zenith angle. Since Bye, and Lsky are both 


constant, the path representing the average sea radiance is the same path as that for the 


average reflectance. 
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From the Fresnel equations, the E-vector amplitude reflection ratios are : 


_ n, cos, —n, cos6, 
ioe a (C.4a) 
n,cos@, +n, cos8, 


_ n, CoS, —n, cosG, (C.4b) 


5 
n,cos@, +n, cos6, 


Taking n;=1 and n, =1.33 and making use of the Snell's law, we can find the reflectance, 


p as a function of 8; , where : 
(0; )=n? +47 (C.5) 


Hence, the average reflectance is : 


7 2 (C.6) 
(p) == Jo(@)46 


Again, performing the integral numerically, we found the average reflectance of the water 
surface to be 0.208. This value corresponds to the slant path to sea with a zenith angle of 


118.1 degree. 
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APPENDIX D. LIST OF PROGRAMS 


File Name Function 


procedure | Main procedure to extract meteorological data and 
determine aspect angle and pixel size 


Reads the local time (PST) from *.pte files 


Reads meteorological data and returns a met 
vector containing all the required information 


Start.pro 


function 
function 


Readtime 


Readmet.pro 


function 


Planes.pro function Selects the planes that are visible for a given 
aspect angle 


Leia_pte.pro Extracts the images, information or time vector 


procedure | Calculates the geometry inputs for SeaRad files for 
each pixels : 


Writes SeaRad input files between two altitude 


Planein.pro 


Seain.pro procedure 


Writes SeaRad input files for average sea and sky 
paths 


Avein.pro procedure 


Posn.pro ction Calculates the position of each node 


Rposn.pro function | Calculate the position of each node after rotation 
to a given aspect angle 


Vector.pro function Calculates the vector, V, from the sensor to the 
node position 

Vsea.pro function Calculates the vector from the nodes to the sea 
Calculates the magnitude of the V vector 


Vmag.pro ction 


| 


Zenith.pro function Calculates the zenith angle at the sensor for each 
node 


Part2 
Seal 


i 
© 
7) 


a 
r 

© 

Nn 


Copy all input files in Case1/Inputs directory to 
batch file | tape 5 and call X1 


Execute SeaRad and copy output file to 
Casel/Outputs directory 


batch file 










Continue.pro procedure | Main procedure to continue the analysis after 
running SeaRad 
Readtemp.pro | function Reads the temperature matrix 


Seaout.pro function | Reads the output files from SeaRad and returns a 
matrix with the transmittance and path radiance 

Aveout.pro function | Reads the output files from SeaRad and returns a 
vector of average sky and sea radiance 


function Calculates the in-band radiance 


ction | Aligns all the radiance values of planes into one 
matrix 


Block.pro procedure | Determines areas of overlap between planes 
Overlap.pro ction Corrects the radiance values in the overlapped area 


Radtot.pro function | Calculates the equivalent blackbody temperature 
from radiance 


Convtemp.pro | function Converts radiance matrix to temperature 


Ttotl.pro function Calculates thermal level from temperature based 
| on calibration constants 
Ditotl.pro function Converts digital level to thermal level 
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APPENDIX E. STRUCTURE OF CONVERTED DATA FILE 


The converted files have the following properties : 

- all polarisation cases for one test case are combined into one file rather than 
having two separate files | 

-each image is a 125 x 256 integer array 


- all essential information is put into a single integer matrix 
The structure of the file and the format for header information is tabulated below. 


Table E1. Structure Of A Basic (.Pte) File 


Name of the original (.PTW) file containing string(20) — 3 
horizontal and vertical polarisation cases | 












Table E2. Information Vector 


















«| Mint Gt a se 
5 


an a 
5 [Seconds tame) 







Hundreds (first frame) 


8 Ramtec ofpietsprtne 
9 [Namberofins dC 


Number of frames : horizontal polarisation 
Number of frames : horizontal polarisation 
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APPENDIX F. LOCATION OF MEASUREMENTS 
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APPENDIX G. STORAGE FORMAT FOR METOC DATA 


There were two meteorological towers installed onboard R/V Point Sur for 
recording the meteorological data. The storage format for the two towers are slightly 
different as the METOC 2 system has two additional temperature and relative humidity 


sensors. The two formats shall be described below. 
A. METOC1 DATA 


These data files are identified by the filename "eop*.txt" where * is the GMT day 
of the year. Each file is started at 00:00 GMT and the data are averaged every 30 seconds 
and recorded. For each line of a record, the various readings from the sensors are 
separated by spaces across a line of text. At the end of the record, a carriage return and 

line feed character starts the next record on a new line. The order of the sensor readings is 


as follows : 


01 number of readings in this average, nn 

02 year, yy 

03 Julian date, ddd 

04 hours &minutes, hhmm 

05 seconds, ss 

06 relative wind speed (m/s), mmm 

07 relative wind direction (degrees), ddd 

08 ship speed (knots), kk.k - information not available 
09 ship direction (degrees), ddd.d - information not available 
10 true wind speed (m/s), mmm 

11 true wind direction (degrees), ddd 

12 Air temperature (degrees C), tt.t 

13 relative humidity (%), pp 


14 pressure (millibars), pppp.p 
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15 sea surface temperature (degrees C), tt.t 

16 GPS time (GMT), hhmmss 

17 GPS latitude (degrees & minutes), ddmm.m 

18 GPS north/south indicator, either 'N' or 'S' 

19 GPS longitude (degrees & minutes), , ddmm.m 

20 GPS east/west indicator, either 'E' or 'W' 

21 GPS speed over ground (knots), kk.k - information not available 

22 GPS course over ground (degrees), ddd.d - information not available 


23 GPS antenna height (m), mm.m 


A sample of the data is below : 


30 1996 99 2352 32 2 138 O 090 5 288 16.5 — 75 1016.9 
16.8 000035 3311.2 N 11728.2W 6.3 124.3 46.3 


B. METOC 2 DATA 


METOC 2 data files are identified by the filename "eop30sec.* " where * is the 
GMT day of the year. The recording format is the same as METOC 1 files except that 


there are 4 additional sensor readings : 


Q1 number of readings in this average, nn 

02 year, yy 

03 Julian date, ddd 

04 hours & minutes, hhmm 

05 seconds, ss 

06 relative wind speed (m/s), mmm 

07 relative wind direction (degrees), ddd 

08 ship speed (knots), kk.k - information not available 

09 ship direction (degrees),.ddd.d_ - information not available 
10 true wind speed (m/s), mmm 


11 true wind direction (degrees), ddd 
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12 Air temperature (degrees C), tt.t 
_ 13 relative humidity (%), pp 
14 pressure (millibars), pppp.p 
15 sea surface temperature (degrees C), tt.t 
16 Air temperature A (degrees C), tt.t 
17 relative humidity A (“%), pp 
18 Air temperature B (degrees C), tt.t 
19 relative humidity B (%), pp - 
20 "***" | Campbell-GPS data delimiter 
21 GPS time (GMT), hhmmss 
22 GPS latitude (degrees & minutes), ddémm.m 
23 GPS north/south indicator, either 'N' or 'S' 
24 GPS longitude (degrees & minutes), , ddmm.m 
25 GPS east/west indicator, either 'E’ or "W' 
26 GPS speed over ground (knots), kk.k - information not available 
27 GPS course over ground (degrees), ddd.d - information not available 


28 GPS antenna height (m), mm.m 


A sample is below : 


30 1996 99 2352 32 2 138 0 090 5 288 165 75 1016.9 
16.8 0.0 0.0 0.0 0.0 *** 000035 3311.2N 11728.2W 6.3 124.3 


46.3 
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Fileno. | Day of | GMT | Aspect | Range | Seatemp. | AMP | Current | Average Wind Line of 
(base*.pte) | year angle (km) (°C) wind wind azimuth Sight 
(deg) speed speed relative to (LOS) 
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APPENDIX I. NAVY AEROSOL MODEL AND AIR MASS 
PARAMETER 


A. NAVY AEROSOL MODEL [Ref. 5] 


The US Navy Aerosol Model (NAM) was developed to compute the aerosol 
extinction coefficients for visual and infrared wavelengths over a marine boundary 
environment. It is built from meteorological and aerosol size spectra databases obtained 


near the surface layer of several oceanic environments. 


The aerosol size distribution of the NAM is attributed to 3 distinct aerosol 
production mechanisms. The first mechanism is air mass driven while the second and 
third mechanisms are wind driven. Correspondingly, the NAM aerosol size distribution is 


a sum of 3 independent log-normal distribution functions. This is given as follows : 


nye SE sent fol El) an 


where n(r) is the number of particles per cm? per radius r, in microns (um), 
k is the visibility scaling factor, 
r, are the modal radii ( 0.03pm, 0.24 um and 2.0um) measured at 80% RH, and 
Aj = 2000 P’ , where P is the air mass parameter, 
A> =5.866(Vave - 2.2) , where Vavye is the 24-hour average wind velocity, 


A3; =10 ae , where V is the instantaneous wind velocity, 


f = [(2-Rh/100)/6(1-Rh/100) 7°" 
B. AIR MASS PARAMETER [Ref. 5] 


Of the 3 amplitudes A, to A3 given above, the largest and the most difficult to 
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measure is the first one, that is, the air mass driven component. The air mass parameter, P 
is an integer ranging from 1 for clean oceanic air to 10 for polluted continental or coastal 


air. 


There are several methods to obtain the air mass parameter. These methods range 
from measuring the atmospheric Radon-222 gas concentration to computing air mass 
trajectories to remote sensing using satellite. Each method has its advantages and 


disadvantages. 


Traditionally, the Radon-222 gas concentration method is the most reliable. 
Radon-222 is an inert radioactive gas with a half-life of 3.82 days produced by the decay 
of Radium-236. Radium-236 is found in various quantities in the soils and rocks of the 
earth's crust. After formation in the ground, the radon diffuses into the atmosphere and 
becomes attached to aerosol particles that are transported by the prevailing wind over the 
ocean. Since radon is mainly produced on land, it is used as an air mass tracer to 
characterise the air mass history. The air mass parameter, P is then simply given by the 


equation : 
P = INT =| +] (1.2) 


where INT truncates to the integer value and Rn is the atmospheric Radon-222 gas 
concentration in picocurries per cm’. Strictly speaking, the concentration of the radon in 
different areas will not be the same. This is because it depends on the soil conditions, 
atmospheric conditions, radium deposit in the crust and strength of the radium source. 
Therefore, it is the relative change in radon gas concentration that is meaningful. 
Equation I.2, however, uses the absolute radon concentration. This is because the factor 
of 4 was established for the region at San Nicolas Island. This factor will have to be 


modified for other geographical regions. 
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APPENDIX J. SEARAD INPUT 


SeaRad input requires the submission of a series of a minimum of 5 cards. The 
number of cards required depends on the mode of execution. There are a total of 19 cards 
of which 5 are the main cards and the rest are optional cards. A complete description of 
each card can be found in the file "Modtran10" that comes with program. However, this 
does not include the modifications done in SeaRad. The file "notes.txt" logs all the 
modifications and should be read in conjunction with the original card description. This 
appendix will only explain the parameters that are required in the cards needed in this 


thesis. 


CARD 1 (line one) - Format (1315, F8.3, F7.2) 


MODEL, ITYPE, IEMSCT, IMULT, M1, M2, M3, M4, M5, M6, MDEF, IM, NOPRT, 
TBOUND, SALB 


F - use LOWTRAN 7 
MODEL - selects one of the 6 geographical model atmospheres or specifies that 
user-defined meteorological data are to be used 

: mid-latitude summer (2) was used 

ITYPE - indicates the type of atmospheric path 
: slant path to space (3) was used for average sea and sky radiance 
calculations and slant path between two altitudes (2) was used for pixel- 
sensor atmospheric compensation 

IEMSCT - selects the execution mode of the program 
: radiance mode with solar scattering (2) used 

IMULT - determines execution with multiple scattering 


: multiple scattering (1) was used 
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M1 ,.M2,M3,M4,M5,M6 - parameters to modify or supplement the altitude 
profiles of temperature and pressure, water vapor, 


ozone, CH4, N2O,CO 


> default was used (0) 
MDEF - use default profile for C02, 02, NO, S02, N02, NH3, HNO; 
: default was used(0) 
IM - additional switch for user defined atmosphere 


: set to normal operation (0) 
NORPT - switch to minimize printing of transmittance, radiance and atmospheric 
profile tables 
: set to minimise (1) 
TBOUND - boundary temperature in K 
: set to sea temperature. Note : the sea temperature is only used if the path 


intercepts the earth 


SALB -. surface albedo of the earth at the location and average frequency of the 
| calculation (0.0 to 1.0) 
: set to 0 
CARD 2 (line2) —_- Format (615, 5F10.3) 


IHAZE, ISEASN, IVULCN, ISCSTL, ICLD, IVSA, VIS, WSS, WHH, RAINRT, 
GNDALT 


[HAZE - selects the aerosol model for extinction calculation and default 
meteorological range 
: Navy aerosol model used (3) 
ISEASN - selects the seasonal dependence of tropospheric and stratospheric profile 
: set to default (0) i.e. model selected in card I 
IVULCN - selects the aerosol profile above the stratosphere to 100 km 
: default stratospheric background used (0) | 
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ICSTL - air mass character used only with Navy aerosol model (1 for open ocean 
to 10 for strong continental influence) 
: calculated values from radon gas concentration used (see Appendix F) 
ICLD - specifies cloud model and rain rates 


: no clouds or rain (0) used 


IVSA - determines the use of Army vertical structure algorithm for aerosols in 
boundary layer 
: not used (0) 
VIS - specifies meteorological range 
: default used (0) 
WSS - specifies current wind speed in m/s 


: meteorological data used (see Appendix F') 
WHH - 24 hour average wind speed in m/s 
: meteorological data used 


RAINRT - rain rate in mm/hbr 


: default used (0) 

GNDALT - altitude of surface relative to sea level in km 
: setto 0 

CARD 3 (line 3) - Format (6F 10.3, [5) 


H1, H2, ANGLE, RANGE, BETA, RO, LEN 


HI - initial altitude in km (in radiance mode this always defines the altitude of 
the sensor) 
: 0.017 km for pixel-sensor path and pixel height for sea and sky path 
H2 - final altitude in km 
: pixel height for pixel-sensor path and 0 for sea and sky path 
ANGLE - zenith angle measured from H1 in degrees (between 0 to 180), 0 degrees 


is looking straight up 


71 








: calculated for each pixel 
RANGE - path length in km 


: not needed as the first parameter is sufficient to specify the geometry 


BETA - earth center angle subtended by H1 and H2 in degrees 

: not needed as the first parameter is sufficient to specify the geometry 
RO - radius of the earth at the particular geographical location in km 

: default radius of earth i.e. 6371.23 km (0) used 
LEN - select normal or long path through tangent height calculations 


: normal operation (0) used 
PSI - azimuth of wind direction in degrees measured positive East of North 
from observer line of sight 
: meteorological data used 
Seaswitch - switch for turning sea modification on or off 


: set to True (T) 


CARD 3A] (line 4) - Format (8F 10.3) 
(required when IEMSCT =2 ) 


IPARM, IPH, IDAY, ISOURC 


IPARM - controls the method of specifying solar/lunar geometry on CARD 3A2 
: observer lat. and long. used (1) 

IPH - determines the type of phase function to be used 
: MIE generated database for LOWTRAN models (2) used 

IDAY - day of the year in GMT | 


: calculated from meteorological data 
ISOURC - type of extra-terrestrial source to be considered 


: sun (0) selected 
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CARD 3A2 (line 5) - Format (8F10.3) 


PARMI1, PARM2, PARM3, PARM4, TIME, PSIPO, ANGLEM, G 


PARMI1 - observer latitude in degrees (-90 to +90; + means North) 
- 32.917 degrees | 

PARM2 - observer longitude in degrees (0 to 360) 
: 117.517 degrees 

PARM3 - not used 

PARM4 - not used 

TIME - GMT time in decimal hours 


: calculated from meteorological data 
PSIPO - path azimuth in degrees East of North 
: calculated from meteorological data 


ANGLEM _ - phase angle of the moon 


: not required 

G | - asymmetry factor for use with HENYEY-GREENSTEIN phase function 
: not required 

CARD 4 (line 6) - Format (3F10.3) 


IV1, 1V2, IDV, IFWHM, IFILTER 


IVi - initial frequency in wave number in cm” 
: 830 cm" i.e. 12 um used 
[V2 - final wave number in cm” 


: 1250. cm” i.e. 8 ym used 
1 


IDV - frequency increment in cm 
- 10 cm" used 
IFWHM - frequency band at full width half maximum 
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: 5 em” used 
IFILTER - specifies the filter to be applied 
: no filter (0) used 


CARD 5 (line 7) - Format (15) 

IRPT 

IRPT - determines what program will do next 
: set to end program (0) 


A sample file is shown below. 


F2 3 2 1 0 0 0 0 90 


0 


3 0 0 5 O 9 0.000 7.000 
0.017 Q.001 90.531 0.000 0.000 


0 


1 2100 0 

32.660 117.240 0.000 0.000 21.181 
830 ©1250 10 5 

0 


74 


0 O 
5.200 
0.000 


245.016 


1 287.95 
0.000 . 
0 50.984 


0.000 


0.00 
0.000 
T 


0.000 





APPENDIX K. SEARAD OUTPUT 


The format of the SeaRad output files depends on the type of calculation that is 
done. Although all cases were executed in thermal plus solar radiance mode, their output 
is slightly different because of the way the path geometry was specified. A sample output 


of each case 1s below. 


A. OUTPUT FOR AVERAGE SEA RADIANCE PATH 


*#*#* SEARAD, A MODIFICATION OF LOWTRAN7 ***** 
DATE: 02/08/1999 | TIME: 15:09:32.87 
THERMAL PLUS SOLAR RADIANCE MODE 
MULTIPLE SCATTERING USED 
MARINE AEROSOL MODEL USED 
WIND SPEED 6.00 M/SEC 
WIND SPEED 3.75 M/SEC, 24 HR AVERAGE 


RELATIVE HUMIDITY = 76.11 PERCENT 
AIRMASS CHARACTER = 2.0 


VISIBILITY = .00KM 
SLANT PATH TO SPACE 
Hl =  .011 KM 
HMIN = .000 KM 
ANGLE= 118.1001 DEG 
FREQUENCY RANGE 
[V1 = 830 CM-1 ( 12.05 MICROMETERS) 


[V2 = 1250CM-1 ( 8.00 MICROMETERS) 
IDV = 10 CM-1 

IFWHM = 5 CM-1 

FILTER = v 
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SUMMARY OF THE GEOMETRY CALCULATION 


011 KM 

H2 =  .000 KM 
ANGLE = 118.100 DEG 
RANGE = 0.023 KM 


= 
| 


BETA =  .000DEG 
PHI = 61.900 DEG 
HMIN =  .000 KM 
BENDING= = .000 DEG 
LEN = 0 


SEA AT 287.45 K REPLACES BLACK BODY BOUNDARY 
UPWIND = -5.006 DEG EAST OF LINE OF SIGHT 


SUMMARY OF OBSERVATION GEOMETRY 


BETA = .00018 DEG 

PATH AZIMUTH = 314.994 DEG EAST OF NORTH 

WIND AZIMUTH = 309.988 DEG EAST OF NORTH 

RECEIVER LATITUDE = 32.917 NORTH OF EQUATOR 
RECEIVER LONGITUDE = 117.517 WEST OF GREENWICH 
FOOTPRINT LATITUDE = 32.917 NORTH OF EQUATOR 
FOOTPRINT LONGITUDE = 117.517 WEST OF GREENWICH 
SUBSOLAR LATITUDE = 7.763 DEG NORTH OF EQUATOR 
SUBSOLAR LONGITUDE = 179.866 DEG WEST OF GREENWICH 


VALUES SEEN FROM FOOTPRINT 


RECEIVER ZENITH ANGLE= 61.900 DEG 


RECEIVER AZIMUTH = 174.994 DEG WEST OF UP WIND 

SOLAR ZENITH ANGLE = 62.650 DEG 

SOLAR AZIMUTH = 48.828 DEG WEST OF UP WIND 

SOLAR SPECULAR TILT = 40.741 DEG ( 7.03 SIGMA, PROB = 1.720E- 


10) 
ZERO RANGE UNPOLARIZED VALUES 


29.68807 W M-2 SR-1 (AV. EMISS. .9426) 


SEA EMISSION = 
SKY REFLECTION = 1.05408 W M-2 SR-] 
SUN GLINT = .00000 W M-2 SR-1 
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TOTAL RADIANCE = 30.74216 W M-2 SR-1 
BLACK BODY TEMP.= —_13.0C 


FULL RANGE UNPOLARIZED VALUES 


SEA EMISSION 29.14501 W M-2 SR-1 

SKY REFLECTION 1.03208 W M-2 SR-1 

SUN GLINT = .00000 W M-2 SR-1 

PATH TO FOOTPRINT = .64978 W M-2 SR-1 (AV. TRANS. .9792) 


| 


TOTAL RADIANCE = 30.82687 W M-2 SR-1 


BLACK BODY TEMP. = 


(WM-2SR-1) (WM-2SR-1) (%) 
ZERO RANGE POLARIZED VALUES 
SEA EMISSION 14.26001 15.42806 -3.9 
SKY REFLECTION 85931 17277 66.5 
SUN GLINT 00000 00000 0 
TOTAL RADIANCE 15.13789 15.60427 -1.5 
BLACK BODY TEMP. (C) _ -22.0 -20.7 
FULL RANGE POLARIZED VALUES 
SEA EMISSION 14.00002 15.14499 3.9 
SKY REFLECTION 3.10762 997732 16.6 
SUN GLINT 00000 00000 0 
PATH TO FOOTPRINT 32489 32489 0 
TOTAL RADIANCE 15.18422 15.64265 15 
BLACK BODY TEMP. (C) _ -21.9 -20.6 


B. OUTPUT FOR PIXEL TO SENSOR PATH 


13.2 °C 


HORIZONTAL VERTICAL 


(H-V)(B+V) 


#eKEE SEARAD, A MODIFICATION OF LOWTRAN7 ***** 


DATE: 02/08/1999 


T7 


TIME: 15:09:10.40 








THERMAL PLUS SOLAR RADIANCE MODE 

MULTIPLE SCATTERING USED 

MARINE AEROSOL MODEL USED 
WIND SPEED = 6.00 M/SEC 
WIND SPEED = 3.75 M/SEC, 24 HR AVERAGE 
RELATIVE HUMIDITY = 76.11 PERCENT 
AIRMASS CHARACTER = 2.0 
VISIBILITY =  .00KM 


SLANT PATH, H1 TO H2 


Hi = .017KM 
H2 = .011 KM 
ANGLE = 90.367 DEG 
RANGE = ~ .000 KM 
BETA =  .000 DEG 
LEN = 0 
FREQUENCY RANGE 
Vi = 830 CM-1 ( 12.05 MICROMETERS) 
IV2 = 1250CM-1 ( 8.00 MICROMETERS) 
IDV = 10 CM-1 
IFWHM = 5 CM-1 
IFILTER = 0 


SUMMARY OF THE GEOMETRY CALCULATION 


HI =  .017KM 

H2 = .O11KM 
ANGLE = 90.363 DEG 
RANGE = ~~ .948 KM 
BETA =  .009 DEG 
PHI = 89.642 DEG 
HMIN =  .011 KM 
BENDING= —= .000 DEG 
LEN = 0 


TBOUND SET TO .10 K FOR MARINE SKY 
INTEGRATED ABSORPTION = 108.01 CM-1 FROM 830 TO 1250 CM-1 
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AVERAGE TRANSMITTANCE = _ .7428 


MAXIMUM RADIANCE — = 3.722E-02 W M-2 SR-1 (CM-1)-1 AT 830.0 CM-1 
MINIMUM RADIANCE = = 1.351E-02 W M-2 SR-1 (CM-1)-1 AT 1090.0 CM-1 
BOUNDARY TEMPERATURE = 10K 

BOUNDARY EMISSIVITY = _ 1.000 

FILTERED RADIANCE = 8.751E+00 W M-2 SR-1 

BLACKBODY TEMPERATURE =  -44.0C 
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APPENDIX L. STORAGE FORMAT FOR SKIN TEMPERATURE 
DATA : 


The skin temperature of the ship was recorded through 14 thermistors attached to 
selected locations around R/V POINT SUR. The temperatures are read every 20 seconds 
and the original files are identified by the filename "flir13 /*Apr" where * is the day. 


The starting time and day of the recording is 00:00 local time (PST) on the first 
day of each month. At top each file, the starting time and day of the first line of record of 
the file is also indicated. For each line of record, there are 17 entries separated by spaces 
across a line of text. At the end of the record, a carriage return and line feed character 
starts the next econ on a new line. The first entry is the local time in seconds from the 
beginning of the first day of the month, and this is followed by the readings of the 16 
- thermistors in ascending order. Readings 1 and 9 are reference readings and should 


always be at 25.0. All temperatures are in degrees Celsius. 
A sample of the data is below : 


april9.1 

Mon,Day,Year= 4 9 1996 

Hour,Min,Sec= 20 8 27 

Volts, istart,istop,nsam,iref 1.49939 1 16.100 1 


763707 =.25.001 3=14.458 = 17.209 13.710 14.943 16.570 13.946 
14.468 25.001 15.562 14.130 31.901 15.310 27.986 15.783 
15.153 


8] 
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APPENDIX N. ORIGINAL FILES AND R/L SETTINGS 


Original filename | Range setting (R) Level setting (L) 
NSCA0912.ptw 
NSCA0913.ptw 
NSCA0915.ptw 
NSCA0916.ptw 
NSCA0918.ptw 
NSCA0919.ptw 
NSCA0927.ptw 
NSCA0928.ptw 
NSCA0930.ptw 
NSCA0931.ptw 
NSCA0933.ptw 
NSCA0934.ptw 
NSCA0936.ptw 
NSCA0937.ptw 
NSCA0938.ptw 
NSCA0939.ptw 
NSCA0940.ptw 
NSCA0941 .ptw 
NSCA0942.ptw 
NSCA0943.ptw 
NSCA0944. ptw 
NSCA0945.ptw 
NSCA0946.ptw 
NSCA0947.ptw 
NSCA0948.ptw 
NSCA0949.ptw 


Filename 


as 


BaseQ1.pte 


amen 
fs 


as 


BaseQ2.pte 


aw 
— 


mm} Bi Py] mi Rl] By, HB; HY] HY] HP] A 


pod 


Base03.pte 


ped 


Base(04.pte 


— 


Base05.pte 


Base06.pte 


Base07.pte 


jt 


Nn 
bond pa — — nt 


Base(08.pte 


Base09.pte — 


Base10.pte 


i. 


Basel 1.pte 


— 
a 


_ Basel2.pte 


— 


Base13.pte 


pow 





Filename 


Base14.pte 


Basel5.pte 


Basel 6.pte 


Basel 7.pte 


Base27.pte 


Base28.pte 
Base29.pte 
Base30.pte 


Base3 1.pte 


Base32.pte 
Base33.pte 
: Base34.pte 


Base35.pte 


Base36.pte 


Original filename 
NSCA0950.ptw 
NSCA0951 .ptw 
NSCA0952.ptw 
NSCA0953.ptw 
NSCA0954.ptw 
NSCA0955.ptw 
NSCA0956ptw 
NSCA0957.ptw 

_ NSCA0978.ptw 
NSCA0979ptw 
NSCA0980.ptw 
NSCA098 1 .ptw 
NSCA0982.ptw 
NSCA0983.ptw 
NSCA0988.ptw 
NSCA0989.ptw 
NSCA0990.ptw 
NSCA0991 .ptw 
NSCA0993.ptw _ 
NSCA0994.ptw 
NSCA0996.ptw 
NSCA0997.ptw 
NSCA0999.ptw 
NSCA1000.ptw 
NSCA902.ptw 
NSCA903.ptw 
NSCA905.ptw 
NSCA906.ptw 


4) 


oO 


— i pescmecds 
as 
ve) 


——s —d 
ay 
7) 


ON 
poh 
; Cc Ns 
aN 
oe) 


Level setting (L) 
0 
5 


47 
44 
46 
43 


iN 
ON 


© 
as 
IS 


— pwnd 
=~ | 
—“J] 


om) 
aN 
ee) 
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— 
© 
— sn owt — 


\O 


© 
Ln 


on) 
a 
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a) © 
A ae i il i nl OA 
WtEnt 


aS 
\O 


aN 
\O 


>) 
ao 
oN 


ome) 
>; 
Wt oO 


WG 
© 














Original filename | Range setting (R) Level setting (L) 


Base37.pte NSCA911.ptw 
Base38.pte NSCA914.ptw 
Base39.pte NSCA917.ptw 
Base40.pte NSCA920.ptw 

NSCA921.ptw 


Base41 .pte NSCA923.ptw 
Base43.pte NSCA929.ptw 
_ Base44.pte NSCA932.ptw | 
. NSCA933.ptw 


Base45.pte NSCA935.ptw — 
- Base49.pte | NSCA947.ptw 
_ Base50.pte NSCA958.ptw 
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APPENDIX O. AN EXTRACT ON PERFORMING COMPLEX 
THERMAL MEASUREMENTS FROM THE AGA MANDAL [Ref.9] 


Section 10 of the AGA operating manual describes the thermal measurement 
techniques for making measurements with the AGA Thermovision infrared scanner. The 
extract in the next three pages describes the theory and how complex thermal 


measurements are made with the AGA scanner. 
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SECTION 10 —- THERMAL MEASUREMENT TECHNIQUES 


area) is just Orightened up by the isotherm. For 
ease of operation this is often carried out with the 
PICTURE MODE switch set to INVERTED or BLACK. 


{c) Adjust the second isotherm (LEVEL 2) to indicate 
the colder of the two areas of interest. If this cannot 
be achieved. increase the THERMAL RANGE one 
step at a time and repeat (a) and (b) and {c) again 
until both colder and warmer areas are covered by 
the same selected range. 


NOTE: For overall accuracy maintain the THERMAL 
RANGE sefting at its lowest level consistent with 
displaying both the reference and abject area 
isotherms. 


Without readjusting THERMAL RANGE or LEVEL 
orcceed as foilows: 


iG) Carefully adjust the two isotherm controls ‘o 
indicate the object area to be measured and the 
reference area respectively. This is shown by the 
areas being displayed as saturated wnhite. 


‘e) Read the isotherm scale values at the top edge of 
both isotherm markers, calculate the uifference 
and muitipiy by the THERMAL RANGE setting isee 
figure 2). 


This results in a thermal difference value _.i,, tnat can 
be applied directly to the calibration curve since we 
have postulated ¢ 5=¢€,=1 and 7 g=1. Proceed as 
follows: 


if) Mark the reference temperature on the calibration 
cnart. Read via the curve the corresponding 
thermal value i,. Add the thermai difference vaiue 
—ig- Obtained in (e) to |,. This gives the blackbody 
thermai value |, of the cbject, which is then 
converted to temperature t, by means of the curve. 


NOTE: {f the reference should have higher temperature 
than the object, the thermal difference lop = -!r 
becomes negative and should consequently 3e 
subtracted from the reference thermal vaiue i,. The 
iower half of the isotherm scale has negative readings. 


itis therefore important to note the correct sign of the 


isotherm scale readings when calculating lor 


COMPLEX THERMAL MEASUREMENTS 


The previous instructions have assumed the ideal 
situation, that the object is a perfect blackbocy and tnat 
no external factors influence the measurements. 


in compiex thermal measurement situations, however, 
where several factors influence the measurements, ne 
true object temperature has to be derived dy 
calculation. 


An expression 's set up to include the -adiation that 
reaches the scanner from the opoject. from this 
expression the unknown parameter is 2xtracted ‘e.g. 
abiect temperature or abject emissivity). Before cnis 
derivation is carried out the following variabies shouid 
3e.considered. 


Variables affecting measurement 


The most important external influencing factors are: 


AGA Thermavision 780 
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Emissivity - Actual objects are seldom “black” Tre 
emissivity factor must therefore often be taken into 
account in infrared temperature measurements 35 
detailed in the measurement formulae. !ndivicuai 
object emissivity can be measured or can be ‘ound in 
tables, (e.g. Table 1 page 11.8). 


Ordinarily, object materials and surface treatments 
exhibit emissivities ranging from approximately 0.1 to 
0.95. A highly polished (mirror) surface falls below 0.1, 
while an oxidised or painted surface has greatly 
increased emissivity. Oil-based paint, regardless of 
colour in the visible spectrum, has an emissivity over 
0.9 in the infrared. Human skin exhibits an emissivity 
close to 1 which is used of in medical thermography for 
direct body-temperature readings. 


Surroundings — The reflectivity factor of an cpaaue. 
diffuse surface is ¢ =1- €. 


A low emissivity factor means therefore not only :nat 
the emission from the object is lower than ‘rom 3 
blackbody of the same temperature. It also means that 
undesired radiation from the surroundings is reflected 
into the scanner. The measurement formulae °s 
corrected: for reflected radiation by means of 3 zerm 
inctuding I3, which is related to ambient temperature. it 
is therefore important that the object surroundings 
have a homogenous {ambient} temperature and do not 
include hot areas so positioned that their radiation can 
be reflected by the object. Sometimes an 2fficient 
ambient temoerature has to be estimated to take into 
account radiation sources that cannot be removed or 
shielded. 

Object opacity - Objects may 5e more or 'ess 
transparent to infrared radiation e.g. glass anc plastics. 
The radiation recetved from such an object wouia 
include irrelevant transmitted radiation from the 
Oackground. 3esices the efficient emissivity wouid be 
iow and hard to estimate. However, "he spectral 
cnaracteristics of the scanner may be adapted to tne 
object by means of spectrai filtering such that oniv 
spectrum regions where the emissivity is nigh and the 
transmission is zero are utilised. Where this is cossitie 
to do, thermal measurements can be carried out. 
However, the resuitant emissivity anc reflectivity of tne 
object must still oe taken into account if the spectra! 
filtering does not resuit in an effective blackbody. 


Object size — The detector subdtends a certain solid 
angie in the cbject space. !f the object does not cover 
this angle the detector will receive radiation ‘rom tne 
object background. This will tena to iower che 
temperature difference indicatea between the cbiec? 
and Dackground. !t is therefore recommendec ‘nat the 
size of the :mage area, wnere temperature :s to 5e 
Measured, snouid not de less than 2.5 x 2.5 mm as seen 
on the Black/White monitor chassis display screen. 


The Atmosphere - Certain constituents of the atmos- 
phere absorb infrared radiation in the spectral bands 
being used. The most important gases are water 
vapour (H,Q) and carbon dioxide (CO.;. This acsorption 
wil attenuate che infrared radiation, from the ooyect fo 
the instrument. In some cases ‘e.g. not combustion 
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SECTION 10 ~ THERMAL MEASUREMENT TECHNIQUES 


flames) the self-emission of atmospheric gases Is verv 
high and has to be taken into consideration. 


The atmospheric attenuation can be calculated and 
compensated with great accuracy using either 
computer programs or an apprnximative formula. iSee 
Section 10 page 10.6.) 


Measurement formulae 


In general a number of different parameters are 
involved in temperature determinations using 
radiometric methods. Here are listed all parameters 
used in the following description. They are divided into 
three groups of which the !ast is used only for formula 
derivation. 


Radiation source and atmosphere related parameters 
c object temperature ( C} 


temperature of object surroundings 
(ambient) (°C) 


t - reference temperature ( C} 
tam = atmosphere temperature (“C) 


emissivity factor,Os¢€ <1 
(€ =1 fora “blackbody”’} 


= object emissivity 
emissivity of reference surface 


emissivity of object surroundings 
é ,:s considered =i 


7 = atmospheric correction factor for the 
atmosohere between scanner and opject, 
te. atmospheric transmission corrected 
for systems spectral response 


tT = as above for the atmosphere Derween 
scanner and reference 

do = distance scanner to object 1m) 

Ud, = as above to reference !m) 


calibrated thermal value. Refers to ‘he 
instruments calibration function 
'grapnical curve or calculator program) 
which describes the reiation deneen 
object blackbody temperature and instru- 
ment numerical output on snort distance. 
f=fft) ¢tU) 

le = calibrated thermal value fort=t, U) 
calibrated thermal value for t=t, IU) 
calibrated thermal vaiue fort=t, ‘IU) 


_ amas 7 $ — ‘ 
latm = calibrated thermal value fort=tarm I) 


Instrument reading parameters 


i measured object thermal vaiue. Instru-— 
ment reading of radiation received from 
the object surface via the atmospnere i!U) 
10.4 
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l, = as above for the reference (IU) 


The instrument reading consists of twoterms: i’ =L-1, 

where | 

L = THERMAL LEVEL setting of the :nstru- 

ment (IU) 

Lo = as above at object measurement (IU) 

L = as aboveeference measurement (IU) 

= measured relative thermai value = 
instrument scale or display reading x 
THERMAL RANGE setting (IU) 


lo = as above referring ta object measurement 
(WU) 


F = as above referring to reference measure- 
ment (IU) 


lor = as above referring to the difference 
between object and reference measure- 
ment (IU) 


—loron'r 


Radiation notations 


S = received and detected radiation from a 
blackbody of temperature t, at snort 
distance :Photons/sec} 

So = as above when t=t, iPhotons/sec) 

Sy = as above when t=t, i‘Photons/sec! 

Sstm = as above WHEN t=tarm, (Photons/sec} 

So = radiation received via the atmospnere 


fram the object surtace iPhotons/sec! 





ATMOSPMERE 


SURROUNDINGS 


Fig. 3 Radiation conditions in the general 
measurement situation 


Figure 3 iilustrates scnematicaily the radiation 
conditions in the general measurement situation. Note 
that the reflection factor cf the object surface is 1- € 5 
‘or an Opaque object, and that the emission factor of the 
atmospnere is {1- %). With the support of this figure 
the received radiation can be written: 

So = Ta: €g SQ 7 Tgi1 -€5)Sg 511-7 i Sam 
=object radiation+reflected radiaticn- 

atmosphere radiation 


AGA Thermovision 780 








SECTION 10 - THERMAL MEASUREMENT TECHNIQUES 


This radiation relationship can be converted into a 
thermal value relationship by utilising the fact that AGA 
Thermovision” 780 has a linear ,....ton counting 
detector. 


We can therefore write the thermal vaiue 
l=CxsS 
where C is an empirical instrument factor 


; La le: 4 
Replacing S,',S 5. Sa and Sarm by = e a and = 
respectively results in 

yo ig Soo lg, F Tae ee Pa 


This is the general measurement formula describing 
how the measured thermai value |, relates to the three 
radiation sources: object, surroundings and 
atmosphere, together with their associated parameters 
emissivity and transmission. The object is considered 
opaque [1.e. no transmitted radiation througn the 
object) and all the relevant surrounding surfaces are 
assumed to have the same (ambient) temperature f3 
(from that assumption follows automaticaily €,=1) 


The radiation terms are expressed as blackboay 
thermal values determined by the calibration function 
at respective source temperature. 


This formula is a suitable starting point for derivation of 
special measurement formulae adaoted to special 
measurement situations. 


For direct measurements (without reference) the 
formula can be rewritten, using the relation lo =Lo- 


lo: 
Soive for |): 
pe oe fe-— lg att = Mage :2) 
ao fg 69 “9% . 


In case the object temperature ts known out the coject 
emissivity unknown, che latter can be caiculated. Solve 
for ¢€ 9: 


= — « 4 


nat 





ror relative measurements with a xnown temperature 
reference, corresponding formuige are derived as 
follows: 


i = T € | + —_ j + if - } 

‘go “a 'o Te (I €.,/ _ il T,! ae 
i; = Lee —€é)} > {)-T 
t c, ne 0 (1 ©, / Uae ul 7} eee 


were tag and is, refer to the surroundings of abject 
and reference respectively. 


12 =| Pes | 


AS previously described under relative measurements, 
the THERMAL LEVEL shouid be kept constant wnen 
using this method. Thus L, =L, and from that follows: 


AGA Thermovision 780 
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in a practical measurement situaticn, we assume that 
odject and reference are olaced close together. Hence 
TF and bao =lar=l- With these assumotions we 
get 





« 
or e. { = ; : 
ly gates sens + (f!-— 3 (dq) 
0 &o a ES 


Solve for € 4: 


Si ltt = 
€ = ee Re or , rs (5) 
° _ _ 
o's I) 


The above derivated formulae can be a great deal 
simplified under speciai conditions. For instance if all ¢ 
and +r are close to unity we get 


for direct and relative measurement respectively, i.e. as 
assumed in the first part of this section. 


The formula for relative measurement will be reduced 
and thereby the overall accuracy improved, if €,=. €, 
and T)= 7. 


This requirement is met :f for example the reference 
surface 's a part af the object surface itsetf where the 
temperature is known. The formulas then reduced to 


Aj 
[eee heap 
 ] r 
to: Fas, 
ATMOSPHERIC INFLUENCE ON 


INFRARED TEMPERATURE MEASUREMENTS 


The atmosphere :s not fuily transparent ¢o infrared 
radiation. It attenuates radiation passing through it ana 
also emits radiation, tne amcunt dependent on the 
remperature. 


The measurement formuias previously descmbed 
indicated how to correct for the atmospheric influence 
provided that the correction factor r is known. 7 
depends on a number of parameters. The most 
important darameters are iisted Delow together with 
the actual values ‘or what ‘s defined as ‘standard 
atmosphere’. (See Table 1.) 


Table 1 Parameters influencing the 
atmospheric correction factor 7 


a ase ee eT 
Parameter Standard atmosphere! 





Scanner type ‘SW, 3BAR, LW} 

Distance scanner to obrect | 

Air temperature | 

Air pressure | Tatm=i1013 mb 

Relative humidity ‘H.0) | 50% ; 

Carbon dioxide (CO,j | 300 ppm 
| 


~ 
78 - 
+ A 


Carbon monoxide 'CC) | tsopm | 
Dinitragen oxide (N,O) |; 0.2 ppm : 
Ozone (0,) 9.355 ppm. : 
Methane {CH,) 2 pom 
Aerosol visibility ‘particies) 10 km 


i Fo Ee ene nee 
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